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Abstract. We present a solution to a problem suggested by Philippe Biane: We 
prove that a certain Plancherel— type probability distribution on partitions converges, 
as partitions get large, to a new determinantal random point process on the set Z+ 
of nonnegative integers. This can be viewed as an edge limit transition. The limit 
process is determined by a correlation kernel on Z_|_ which is expressed through the 
Hermite polynomials, we call it the discrete Hermite kernel. The proof is based on a 
simple argument which derives convergence of correlation kernels from convergence 
of unbounded self— adjoint difference operators. 

Our approach can also be applied to a number of other probabilistic models. As 
an example, we discuss a bulk limit for one more Plancherel-type model of random 
partitions. 



Introduction 

This work appeared as our attempt to solve a problem posed by Philippe Biane. 
In [Bi2] he considered a model of random partitions arising from decomposition 
of tensor spaces (C N )® n (the Schur-Weyl duality between representations of the 
symmetric group S„ and the unitary group U(N)). The partitions in question 
have at most N nonzero parts which sum to n, and the weight of a partition A is 
proportional to the product of dimensions of the irreducible representations of S n 
and U(N) indexed by A. 

Biane discovered that as n and N go to infinity so that y/n ~ cN then the 
boundary of the Young diagram associated to the random partition A, suitably 
scaled, tends to a nonrandom limit curve given by an explicit formula. The limit 
curve depends on the parameter c > 0. 

If n is fixed while N — > oo then the model turns into the well-known Plancherel 
model of random partitions of n. This agrees with the fact that as c approaches 
0, Biane's limit curve turns into the celebrated Vershik-Kerov-Logan-Shepp limit 
curve for the Plancherel model found in [VK1], [VK2], [LS]. 

Biane's formulas show that the value c = 1 is special: The tangent line to the 
limit curve at one of its endpoints has (in appropriate coordinates) slope —1 for 
c < 1, for c = 1, and +1 for c > 1. Biane's question concerned the local structure 
of the boundary of the random Young diagram at c = 1 near this point of the limit 
shape. 

We address this question in a modified form. Namely, we replace the initial prob- 
ability distribution on partitions by its poissonization with respect to parameter n. 
This procedure is well known, we explain it in §4. One expects that poissonization 
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does not affect the asymptotics but we leave the discussion of this issue out of this 
paper. 

After poissonization, the probability distribution lives on all partitions A = 
(Ai, A2, . . . , Ajv) with at most N nonzero parts and without any constraints on 

|A| = Ai H h Aat. It is convenient to interpret A as an A^-point configuration on 

Z+ = {0,1,2,...} via 

A i-> {xi, . . . ,xn}, Xi = Xi + N — i. 

The weight of A now depends on the poissonization parameter v > (which replaces 
n) and has the form 

const • jj v 1 y • ] [ fa -Xjf. 

i=\ X% ' l<i<j<N 

This is the so-called Charlicr orthogonal polynomial ensemble. 

Our main result is the following statement (see below Theorem 4.1 and Propo- 
sition 3.3). 

Theorem. Fix an arbitrary s e R. Let N = 1,2,... and assume that the parame- 
ter v = v{N) depends on N in such a way that 

v = N 2 + (s + o(l))7V 3 / 2 , N ^00. 

As N — > 00, the probability distribution of{xi, . . . ,xjy} converges to a probability 
measure on 2 z + 7 the set of all point configurations Y on Z + . The correlation 
functions of the limit measure have the form (k = 1, 2, . . . ) 

Prob{F D {t/i, . . . ,y k }} = det[K s (y i ,y j )} 1 < iJ < k , y 1} ...,y k e Z+, 

where, for i,i/eZ + , 



K s (x,y) = {wx\y\2 x +y)-i / e~* H x {t)H y {t)dt 

Js/V2 



/V2 

2 xH^U-%) ■ HJ-%) - HJ-%) ■ yH v _A-%) 

= {-Kx\y\2 x+y )-^ e-^r ^ V ^ ^ - ^ 

x-y 

= UnxWr+vyi e-4 ^ +l( 75 ) ^ ( 7f ) -^ ( 7f ) ^ +l( 75 ) 

x-y 

Here H m is the classical Hermite polynomial, see [KS] . 

The introduction of the additional parameter s above is also due to Biane. 

The determinantal structure of the correlation functions means that the limit 
measure belongs to the class of determinantal random point processes which arise 
in a variety of probabilistic models, see, e.g., [Sol], [So2], [BHKPV], [Ly]. In partic- 
ular, the determinantal processes arise in connection with the Plancherel measure, 
see [Jol], [BOO]. 

We call K s (x,y) (the kernel of the determinantal formula) the discrete Hermite 
kernel. Many similar examples of correlation kernels are known, however, to our 
best knowledge, the discrete Hermite kernel is new. 
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To prove the theorem we have to check that the Charlier correlation kernel (which 
is essentially the iVth Christoffcl-Darboux kernel for the orthogonal Charlier poly- 
nomials) converges to the discrete Hermite kernel. Usually such facts are verified 
using asymptotics of orthogonal polynomials (see, e.g., [Jol] for a different limit 
regime for the Charlier kernel). Such an approach is applicable to our problem 
as well. However, we take another path and extract the needed convergence from 
an abstract theorem concerning strong resolvent convergence of unbounded self- 
adjoint operators. These self-adjoint operators appear as difference operators on 
Z + associated to the Charlier polynomials. This approach seems to be new 1 , and 
it appears to be much less technical comparing to the traditional one. 

To demonstrate the effectiveness of this approach we apply it to another model 
of representation-theoretic origin. This model appeared in the works of Biane [Bil] 
and of Pittel and Romik [PR] , it turns out to be related to the so-called Krawtchouk 
orthogonal polynomial ensemble. In §5 we sketch a proof of the convergence of the 
Krawtchouk kernel to the discrete sine kernel. This result cannot be viewed as new: 
it can be extracted from [IS] , [BKMM] . The point is that our argument is short and 
direct. We show that the result allows one to predict the form of the limit shape 
obtained in [Bil] and [PR]. 

Acknowledgement. We are very grateful to Philippe Biane for posing the prob- 
lem and sharing his insights with us. We are also grateful to Barry Simon for helpful 
advice. The research of A. Borodin was partially supported by the NSF grant DMS- 
0402047. Both authors were also supported by the CRDF grant RUM1-2622-ST-04. 

1. Preliminaries and the Plancherel model 

Let Z + = {0, 1, 2, . . . } be the set of nonnegative integers. Recall that a partition 
is an infinite sequence A = (Ai, A2, . . . ) of nonincreasing numbers from Z + with 
finitely many nonzero terms. The sum of the terms is denoted as |A|. We say that 
A is a partition of n if n = |A|. 

Following [Ma] we identify partitions and Young diagrams. A Young diagram is 
a finite collection of unit squares in the quarter plane with coordinates (i, j), where 
the z-axis is directed downward and the j-axis is directed to the right; a square 
with lower right corner G {1,2,...} x {1,2,...} enters the diagram of A if 

and only if Aj > j. Thus, the diagram A has Ai squares in the zth row with the 
numbering of rows ranging from the top to the bottom, and the total number of 
squares is equal to |A|. The set of all partitions (=Young diagrams) will be denoted 
as Y and the subset of partitions of n e Z + will be denoted as Y„. 

The conjugate partition A' is obtained, in terms of Young diagrams, by trans- 
posing the coordinate axes. Clearly, A^ coincides with the number of nonzero terms 
Ai (or the number of nonvoid rows in the diagram); this number is also denoted as 
£(A). 

The boundary of a diagram A is a broken line going from the point (i, j) — (X[, 0) 
to the point = (0, Ai). It is convenient to add to the boundary those parts of 
the coordinate axes that are below (A' l5 0) and to the right of (0, Ai). The boundary 
of A will be denoted as d\. 

Partitions A e Y can also be regarded as particle configurations on a 1-dimen- 
sional lattice. Let U = Z + 5 denote the set of (proper) half-integers. We assign 



1 Even though the method of deriving the asymptotics of special functions through differential 
equations is well known, we have never seen the same idea applied to correlation kernels. 
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to A an infinite subset of Z': 

X(X) = {\ l - l+ l}, i = 1,2,..., (1.1) 

and we regard X(X) as a configuration of particles sitting at nodes of the lattice Tl . 
The unoccupied nodes of U will be called holes. Note a duality relation between 
particles and holes: reflecting the configuration of holes about we get X(X'). 

The boundary dX (with parts of the coordinate axes included) can be viewed as 
a doubly infinite sequence of horizontal and vertical unit segments, and the particle 
configuration X(X) is a convenient way to encode that sequence. Specifically, a 
node k e Tl is occupied by a particle from X(X) if and only if the the line j — i = k 
meets the boundary dX at the midpoint of a vertical segment. Likewise, the holes 
correspond to the midpoints of horizontal segments. This correspondence makes 
evident the particle/hole duality mentioned above. 

Let S n denote the symmetric group of degree n. The irreducible S^-modules 
are parametrized by the Young diagrams with n squares. Recall the set of such 
diagrams is denoted as Y„. For an arbitrary diagram A € Y„, let dim A denote the 
dimension of the corresponding irreducible module. Equivalently, dim A equals the 
number of standard Young tableaux of the shape A. 

By Burnside's theorem, 

53 (dim A) 2 =n\ 

AGY„ 

The Plancherel measure on Y„, denoted as M ^ lancherel ; i s defined as the probability 
measure with the weights 

M Plancherel (A) = j A G Y„ . 

Regard diagrams A € Y„ as random objects defined on the probability space 
(Y„, M^ lanchclcV ) . As n — > oo, the boundary of the random diagram, scaled by the 
factor of riT 1 / 2 , tends to a (nonrandom) limit curve. This is a well-known result 
due to Logan-Shepp [LS] and Vershik-Kerov [VK1], [VK2] (see also Kerov's book 
[Ke2]). Specifically, if i and j denote the initial row and column coordinates then 
the scaled coordinates are defined as x = i ■ rT 1 / 2 and y = j ■ rT 1 / 2 , and the 
equation of the limit curve in the (x, y) plane can be written as 

x + y = Q(y - x), -2<y-x<2, 

where 

= |(uarcsin| + \J A - u 2 ). 
This result leads to the following conclusions: 

(a) Observe that the limit curve meets the coordinate axes x — and y = 
at points (0, 2) and (2, 0) respectively. This suggests that the first row Ai and the 
the first column X[ of the typical Plancherel diagram A e Y n should grow as 2y / n, 
which is indeed true: see [VK1], [VK2] for a precise statement. Moreover, the same 
holds for each of the largest row and column lengths A^ , X' k , where the index k is 
arbitrary but fixed. 

(b) Let, as above, i and j be the row and column coordinates. Fix u G (—2, 2) 
and let a n be a sequence of positive numbers such that a n — > oo and a n j ' \fn — > 0. 
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Then, as n gets large, the proportion of horizontal (respectively, vertical) steps of 
the random boundary, contained in the strip \j — i — u ■ y/n\ < a n , should be close 
to (1 + fl'(u))/2 (respectively, to (1 — 0'(u))/2), where Sl'(u) is the derivative of 
Cl(u). This statement just means that the slope of the boundary of our random 
Young diagram approximates the slope of the limit curve. 

A somewhat different but essentially equivalent formulation is as follows: Let k n 
be a sequence of half-integers such that k n /y/n — > u e (—2,2). Given n, look at 
the intersection of the line j — i = k n with the boundary of the random diagram 
A e Y n . This is a midpoint of a boundary segment, which can be either horizontal 
or vertical. Then, for n large, the probability to find a horizontal segment should 
be close to (1 + 9!{u))/2. This is indeed true, see [BOO]. 

(c) As u ranges over (—2, 2), the quantity (1 + £l'(«))/2 monotonically increases, 
so that that the probability of finding horizontal fragments increases, too. At the 
cndpoints —2 and 2 the quantity (1 + Q!(u))/2 takes values and 1 (that is, the 
limit curve is tangent to the coordinate axes y = and x = 0). This suggests that, 
typically, each of the differences A fe — A fe+1 or A' fe — A' fe+1 (with k fixed) should grow 
as n — > oo. A much more precise statement can be found in [Ok], [BOO], [Jol]. In 
particular, it turns out that the order of growth of these differences is n 1 / 6 . 

In the next sections we will consider two other models of random Young diagrams 
which may be viewed as deformations of the Plancherel model. 



2. Biane's model 

There is a close relationship between the Plancherel model and the biregular 
representation of the symmetric group. Indeed, the group S n acts on itself by left 
and right shifts. The corresponding representation of the group S n x S n in the 
space of functions on S n has simple spectrum indexed by diagrams A e Y„, and 
(dim A) 2 is just the dimension of the irreducible component indexed by A. Since n\ 
is the dimension of the whole representation space, the Plancherel weight of A is 
equal to the relative dimension of the irreducible component indexed by A. 

Now we apply the same construction but starting with a different representa- 
tion with simple spectrum. Let n and N be two positive integers. Consider the 
tensor space (C^)®" as a bimodulc with respect to the natural commuting ac- 
tions of the symmetric group S n and the unitary group U(N). By the Schur-Weyl 
duality, the representation of the group S n x U(N) in (C") 8 " has simple spec- 
trum which is indexed by Young diagrams A such that |A| = n and £(X) < N. 
Let Y n (N) stand for the set of such diagrams. For A e Y n (N), the dimension of 
the corresponding irreducible component equals dim A • Dimw A, where by Dim at A 
we denote the dimension of the irreducible polynomial U (AT)-module with highest 
weight (Ai, . . . , Ajv). This serves as a prompt for introducing a probability measure 
onY„(A): 

M s r -w cyl(A) = dimA-Dim^A A g ^ (JU) 

(the factor N n in the denominator is the dimension of the whole tensor space). 

Let us take M n c ^ ur_Wcyl as the distribution law for a random ensemble of dia- 
grams A e Y n (N) and ask about the asymptotic properties of this ensemble as n 
and N go to infinity. 

For the first time, this question was addressed by Sergei Kerov [Kel] (see also 
[Ke2, Chapter III, §3]). He showed that if n and N have the same order of growth 
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(that is, n/N tends to a positive constant) then, after scaling by the factor of n~ 1//2 , 
the boundary of the random diagram tends to a limit shape, which is exactly the 
same as in the case of the Plancherel model. This result admits the following 
heuristic explanation: 

If N > n then Y n (N) coincides with Y n , and it is readily checked that 



JV 



lim M^ ur " Wcyl (A) - M n planchml (A), A e Y„ . 



On the other hand, a typical Plancherel diagram A e Y„ has approximately 2y/n 
rows, which explains why the constraint of type £(X) < N — 0(n) turns out to be 
asymptotically negligible. 

Finer results were obtained later by Philippe Biane [Bi2]. 2 He examined a family 
of limit regimes depending on parameter c e (0, +00): 



n — > 00, 



and discovered that, for c fixed, the scaled random diagrams concentrate near a 
limit curve x + y = P c {y — x) depending on c. The curves v = P c (u) are explicitly 
described in [Bi2, §3.1], they provide an interesting deformation of the Plancherel 
curve v — fi(u), which appears as the limit case c = 0. 

Look at the intersection of the curve v = P c {u) with the line v = —it, which 
happens at u = c — 2. A close examination of Biane's formulas (see the end of §3.1 
in [Bi2]) reveals the following fact: 



dP c {u) 



du 



u=c-2 




Our interest is what happens at the critical value c = 1 corresponding to the 
limit regime n ~ iV 2 . 

Assign to A e Y n (N) an iV-particle configuration on Z + : 

X(X) = { Xl ,...,x N }, x l = X t + N-i (2.2) 

(note a difference from (1.1); due to the restriction £(X) < N, A is uniquely deter- 
mined by X(X)). Then the probability space (Y n (N), M^ c ^ ur_Wcyl ) gives rise to an 
ensemble of random TV-particlc configurations on Z + . More generally, we will deal 
with ensembles of infinite random particle configurations as well. Such ensembles 
are examples of what is called a random point process (or random point field). For 
a discrete state space X (in our concrete case X = Z+), a random point process in 
X is determined by specifying a probability measure V on the space Conf(X) = 2 X 
of all subsets in X. Note that Conf (X) is a compact topological space in the natural 
topology. 3 



2 We strongly encourage the reader to look at this paper for a better understanding of what 
follows. 

3 About random point processes in general, see, e.g., [DVJ], [Sol], 
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Conjecture 2.1 (Biane). Consider the random Young diagram A distributed ac- 
cording to the measure M n c ^ ur_Wcyl on Y n (N). Assume that n — > oo and 

N = n 1 ! 2 - \s ■ n 1 / 4 + o(n 1 / 4 ) 

where s is an arbitrary fixed real number. Equivalently, 

n = N 2 + (s + o(l))N 3 / 2 . 

Then the random configuration X(-) converges to a nontrivial random point 
process on Z + depending on s. 

Here convergence means weak convergence of probability measures on the com- 
pact space Conf(Z + ). The limit process is nontrivial in the sense that the limiting 
measure on Conf(Z + ) does not reduce to the delta measure on the empty configu- 
ration or on the configuration coinciding with the whole set Z+. 

In section 3 we introduce the random point processes that appear as limit pro- 
cesses for Biane's model, and in section 4 we verify Biane's conjecture in a modified 
formulation. 

3. The discrete Hermite kernel 

We start with some necessary generalities. Let X be a discrete space and V be a 
random point process in X (that is, a probability measure on Conf(X)). The corre- 
lation functions p n of V (n = 1, 2, . . . ) are probabilities for random configurations 
X e Conf(X) to contain a given finite set {x\, . . . , x n }: 

) = Prob{{xi, . . . , x n } C X}. 

The initial measure V is uniquely determined by the correlation functions p\ , p 2 , ■ • ■ • 
Indeed, for any finite subset A C X, there is a natural projection Conf(X) — > 
Conf (A) given by taking intersection: X i— > Xa = X n A. Let Va be the push- 
forward of V under this projection; this is a probability measure on the finite set 
Conf (A) . Using the inclusion-exclusion principle it is readily seen that Va is deter- 
mined by the values of the correlation functions on A. For instance, for A = {a, b} 
we have 

Prob{JGi = {a,b}} =p2(a,b) 
Prob{Xi = {a}} = pi(a) - p 2 (a, 6), 
Pvob{X A = {b}} = Pl {b) - p 2 (a,b), 
Prob{JGi = 0} = 1 - pi (a) - Pl (b) + p 2 (a, b). 

On the other hand, the initial measure V is clearly determined by collection of the 
measures Va- 

We say that a sequence V\,V 2 , ■ ■ ■ of random point processes in X converges to a 
random point process V in the same space, Vk V, if the corresponding probability 
measures weakly converge. This happens if and only if the correlation functions of 
the processes Vk pointwise converge to the respective correlation functions of the 
process V. Indeed, by the definition of the topology in Conf(X), we have Vk — > V 
if and only (Vk) A — * Va for any finite id, and the latter is clearly equivalent 
to convergence of correlation functions. 
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A random point process in X is said to be determinantal if there exists a function 
K(x, y) on X x X such that the correlation functions are given by the determinantal 
formula 

p n (x 1 , . . . ,x n ) = det[K(xi,Xj)}, n= 1,2,..., 

where the determinant in right-hand side has order n x n. Then K is called the 
correlation kernel of V . Thus, a determinantal process is uniquely determined by its 
correlation kernel. If a kernel K is Hermitean-symmctric, K(x,y) = K(y,x), then 
it serves as a correlation kernel of a random point process if and only if \\K\\ < 1, 
that is, K corresponds to a contractive operator in the Hilbert space £ 2 (X). Indeed, 
this is a very special case of [Sol, Thm. 3]. 

In particular, any projection kernel (that is, the kernel corresponding to a self- 
adjoint projection operator in £ 2 (X)) determines a random point process in X. We 
will introduce now a concrete family of projection kernels which we will need in the 
sequel. 

Consider the semi-infinite Jacobi matrix 

VI 
VT V2 

^Hormitc = \/2 VS . . . (3.1) 

V3 



(the origin of this matrix and of its notation will become clear soon). We agree 
that the rows and columns are indexed by the nonnegative integers x G Z + . Then 
the matrix £) Hermlte determines a symmetric operator in the Hilbert space £ 2 (Z + ): 
by definition, the domain of the operator is the algebraic span of the basis elements 
{S x }, x e Z + . We will denote this operator by the same symbol £) Hermlte . 

Lemma 3.1. The operator £) Hormlto { s essentially self-adjoint. Its closure D H ermite 
has simple purely continuous Lebesgue spectrum. For any Borel set B C R, the 
corresponding spectral projection operator Pb is given by the kernel 

P B (x,y) = (P B S y ,5 x ) = (2TTx\yl2 x +y)- 1 / 2 f e- t2 ' 2 H x (t/V2)H y (t/V2)da, (3.2) 

JB 

where x,y range over Z + and H x (t) denotes the Hermite polynomial of degree x. 

Proof. We will need a few basic facts concerning the classical moment problem: 

Let p be a measure on K with infinite support and with finite moments of all 
orders. Then the space C[t] of polynomials in one variable can be viewed as a 
subspace of the Hilbert space L 2 (R, p). Let C[t] denote the closure of this subspace. 
Finally, let m„ = / t m p(dt) be the moments of p, n = 0, 1, 2, 

(A) The operator of multiplication by t with domain of definition C[t] is essen- 
tially self-adjoint in C[t] if and only if the moment problem associated with the 
sequence {m n } is determinate (that is, p is a unique measure on K with moments 
m n ). See, e.g., [Si, p. 86, Thm. 2]. 

(B) If the moment problem associated with {m n } is determinate then C[t] co- 
incides with the whole space L 2 (R 7 p). See, e.g., [Si, p. 131, Prop. 4.15] or [Ak, 
Cor. 2.3.3]. 



ASYMPTOTICS OF PLANCHEREL-TYPE RANDOM PARTITIONS 



9 



(C) The moment problem associated with {m n } is determinate if the moments 
grow not too fast. For instance, a simple sufficient condition says that the moment 
problem is determinate if the exponential generating series 



ff (z) = ^m„ — 



n=0 



is analytic in a neighborhood of z = (which holds if the function t i— > e zt is 
p-integrable for sufficiently small z). See, e.g., [Si, p. 88, Prop. 1.5]. 

Finally, let Po,pi, . ■ ■ stand for the orthogonal polynomials with respect to p, 
normalized so that J p 2 (i) p(dt) = 1. Then {p n } is an orthonormal basis in C[t]. 
The polynomials p n satisfy a three-term recurrence relation, which means that in 
the basis {p n }, the matrix of the operator of multiplication by t is a (symmetric) 
tridiagonal matrix. 

Now we return to the proof of the lemma. Take as p the normal distribution 

Hermite/^x = _L e^^dt. 



The corresponding polynomials p n are 

H n (t) = (n\2 n )- 1 / 2 H n (t/V2), 

where the H^s are the Hermite polynomials in the standard normalization, sec 
[KS, §1.13]. The three-term recurrence relation for the H n 's has the form 

H n+1 (t) - 2tH n (t) + 2nH n _ 1 (t) = 0, (3.2) 

see [KS, §1.13]. Rewriting this in terms of the HnS we get 

tH n (t) = Vn + lH n+1 (t) + v^ff n -i(t)- 

Thus, the matrix of multiplication by t in the basis {H n } is just the Jacobi matrix 

^Hermite ag defined i n (3.1). 

It is readily checked that pH ™ 1 * satisfies condition (C) above. By (B), the space 
of polynomials is dense in L 2 (R, p Hcrmltc ). Consequently, {H n } is an orthonormal 
basis in L 2 (R, p Hcrmltc ) ; and the correspondence 5 X <-> H x makes it possible to iden- 
tify the Hilbert spaces l 2 (Z+) and i 2 (M, p Hcrmite ). By (A), our operator D Rcrmite is 
essentially self-adjoint. 4 Now it is clear that we have obtained the explicit spectral 
decomposition for the corresponding self-adjoint operator ^Hermite^ ^he closure 
of £) Hermlte . Namely, for any bounded real Borel function \ on K> the operator 
x(-D Hormlto ) is realized as the operator of multiplication by x in £ 2 (R, p Hcrmltc ). In 
particular, taking as \ the characteristic function xb of a Borel set B C R we see 
that the corresponding spectral projection Pg is the operator of multiplication by 
Xb- The matrix of Pb in the basis {H n } is given by formula (3.2), which concludes 
the proof. □ 



4 Of course, all these facts are well known. 
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Definition 3.2. The discrete Hermite kernel with parameter s € R is the above 
spectral kernel corresponding to the set B = [s, oo). That is 

/+oo 
e- t2 / 2 H x (t/V2)H y (t/V2)dt, (3.3) 

where ije Z+. 

The next proposition provides alternative expressions for this kernel: 
Proposition 3.3. For x ^ y, the discrete Hermite kernel can also be written as 

(^UAo*+y\-m P s 2 /2 *g*-i(s/V2) ; H v (s/V2) - H x (s/V2) ■ yHy^js/V^) 

a; - y 

(3.4) 

or equivalently as 

(4 , ,2^-1/2^/2 . H x+1 (s/V2) ■ H y (s/V2) - H x (s/V2) ■ H y+1 (s/V2) 

x-y 

(3.5) 

Proof. The equivalence of (3.4) and (3.5) follows from the three-term relation (3.2). 

Next, we will employ two relations for Hermite polynomials (see [KS, (1.13.6) 
and (1.13.8)]): 

(H n+1 (t))' = 2(n+l)H n (t), (3.6) 
(e-'X(t))' = -e-'X+i(t). (3-7) 

Abbreviating 

C(x,y) = (2irx\y\2 x+v y 1/2 

we have 

r> + 00 



r+oo 

tf s Hermito (x, y) = C(x, y) / e"* ' 2 H x {t/V2)H y {t/V2)dt 

J S 



/■i-oo 

= C(x,y)V2 e-* H x (t)H y (t)dt 

Multiplying by x — y = (a; + 1) — (y + 1) and using (3.6) we get 

(x - y)K?™*%x,y) = ^fi f + °° e-* 2 (H' x+1 (t)H y (t) - H x (t)H' y+1 (t)) dt. 

Integrating by parts and using (3.7) we finally get 

(x-y)K? e ™ ite (x,y) 

= £^) e -^ 2 /2 (h x+1 {8/V2)H v {s/V2) - H x (s/V2)H y+1 (s/V2)) , 

which equals the expression (3.5) multiplied by x — y. □ 
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4. Proof of modified Biane's conjecture 

We will apply a well-known trick called poissonization. Its general idea is to 
make a large parameter n random and obeying the Poisson distribution on Z + with 
large parameter v. Due to the asymptotic concentration of the Poisson distribution 
one believes that the large n limit regime and the large v limit regime are equivalent 
(of course, this claim has to be justified in each concrete situation). On the other 
hand, the latter regime often turns out to be easier to study. 

For instance, as shown in [Jol] and [BOO], application of the poissonization 
procedure to the Plancherel measures M £ lanchorcl leads to determinantal point pro- 
cesses. The same happens for the measures M^ c ^ ur_Wcyl (Lemma 4.2 below). 

By definition, the poissonized version of M^ ur " Woyl , denoted as M P™-Schur-Wcyi ; 
lives on the set 

oo 

Y(N) = (J Y n (N) = {A e Y | t(\) < N} 

n=0 

of all Young diagrams with at most TV rows. This new measure depends on a 
positive parameter v, and is given by 

M Poi SS o„-Schur-Wcyl (A) = g - v ^ M Schur-We y l (A)) A g Y(N). (4.1) 

|A|! 

Clearly, M„ 

oisson c ur oy j g & pj-Q^^ji j^y measure. In the present paper we do not 

justify the poissonization procedure and simply replace the measures M n c ^ ur_Wcyl by 
their poissonized versions in Biane's conjecture. Theorem 4.1 stated below proves 
the conjecture and identifies the limit process. 

Let X Vi n be the random TV-particle configuration on Z + obtained via the cor- 
respondence (2.2) from the random Young diagram A distributed according to the 
measure M P™-Schur-Wcyi Qn y(/V). That is, if {x u ...,x N } = X{\) then 

Prob({x!, . . . , X N }) = M Po^son-Schu r -Wcy 1(A) ^ ^ 

Theorem 4.1. Fix an arbitrary s G M. Let N = 1,2,... and assume that the 
parameter v = v{N) depends on N in such a way that 

v = N 2 + (s + o(l))/V 3 / 2 , TV^oo. (4.3) 

As N — > oo 7 X v /n),n converges to the determinantal point process on Z + with 
the correlation kernel K^ eimltc (x, y) as defined in §3. 

Proof. It suffices to verify that the correlation functions of X V ^ N ^ N converge to the 
respective correlation functions given by correlation kernel K^ CIlnltc (x , y) (see the 
beginning of §3). To do this we will prove that X V ( N ^ N is a determinantal process 
(Lemma 4.2 below) and its correlation kernel pointwise converges to K^ ermltc (x, y) 
(Lemma 4.4 below), which implies the claim of the theorem. 

Consider the weight function for the Charlier polynomials with parameter 9 > 0: 

harlicr (x) = ^ , x € Z+ , 
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see [KS, §1.12]. The iV-particle Charlier ensemble is formed by random iV-particle 
configurations {xi, . . . , Xn} C Z + such that 

N 

Prob({x u ...,x N }) = const(N,6)-l[W^ licr (x l )- J] {xt-xj) 2 , (4.4) 

i=l l<i<j<N 

where const(7V, 9) is the normalization constant (it can be evaluated explicitly but 
we do not need the precise expression). 

Let C m (x; 9) denote the Charlier polynomial of degree m, and let \\C m ( ■ ; 0)\\ be 
its norm in the weighted £ 2 space with the weight function VF ( p harlier : 

oo 

\\c m (. ; e)\\ 2 = J2c 2 m (x-,e)w^(x). 

The normalized functions 

C m (x;6) = (^ charlicr (x)) 1 / 2 ||C m (-;^)ir 1 C m (x;0), m = 0,1,2,... 

form an orthonormal system in the ordinary space £ 2 (Z + ). 

As is well known, the Charlier ensemble is a determinantal point process with 
the correlation kernel 

N-l 

K%y ie *(x,y) - Yl C m (x;9)C m (y;9). 

See, e.g., [Jol], [K6]. This is a projection kernel: the corresponding operator is 
the projection in £ 2 (Z + ) on the iV-dimensional subspace spanned by the functions 
Co, . . . , Civ-i- 

Lemma 4.2. For any v > and N = 1, 2, ... , the random process AVjv coincides 
with the N -particle Charlier ensemble with parameter 9 — v/N . 

Proof. Let us compare the right-hand sides of (4.2) and (4.4). The right-hand side 
of (4.2) is defined by (2.1) and (4.1), this gives 

_„ ^' A ' dim A Dimjv A 

6 |ApnA] • 

We have 



,. , n (xt-xj) 

dim A i<i<j<Af 

~W = n x~\ 

Ki<N 



(4.5) 



(Frobenius' formula, see, e.g., [Ma]) and 

Dimjv A = | J 



7 — I 
l<i<j<N J 



(Weyl's character formula). Finally, 



\X\ = x 1 + --- + x N -^^ 



2 



ASYMPTOTICS OF PLANCHEREL-TYPE RANDOM PARTITIONS 13 



Using these expressions we obtain the right-hand side of (4.4) with an appropriate 
constant. □ 

Given 9 > 0, consider the semi-infinite Jacobi matrix 



Charlicr 






VT 








VI 


-l/Ve 












-2/Ve 










V3 


-3/Ve 



D, 



We also regard it as a symmetric operator in £ 2 (Z + ) whose domain of definition is 
the set of finite linear combination of the basis elements. 



The operator £)C harller - g essen ti a Uy self-adjoint. Its closure Dg haTheT 



Lemma 4.3. 

has purely point spectrum {^7=p | m = 0, 1 
the kernel of the spectral projection operator corresponding to the part of spectrum 



. } . The kernel K^ h g rhcr coincides with 



m 



0,1,...,JV-1 



} 



(4.6) 



Proof. The Charlicr polynomials with parameter 9 satisfy the difference equation 

eC m (x + 1:9)- xC m {x- 9) + xC m (x -l;9) = (9- m)C m (x; 9), 

sec [KS, (1.12.5)]. From the expression for the weight function it follows that for 
the normalized functions C m (x;9), this equation is transformed into 

Vx~+TC m (x + 1;9)- ^f m {x- 9) + Vx~C m (x - 1; 9) = B —^C m (x- 9). 

We see that £)Charher j s precisely the difference operator on Z + standing in the 
left-hand side. 

We claim that £)Charher is essentially self-adjoint and {C m } is the complete set 
of the eigenvectors of the self-adjoint operator £)Chariier i n d ee d, the three-term 
recurrence relation 

-xC m (x; 9) = 9C m+1 (x; 9) - (m + 9)C rn (x; 9) + mC m ^{x- 9) 

(see [KS, (1.12.3)]) and the explicit expression for the norm 

\\C m \\ 2 = e- m e e m\ 

(see [KS, (1.12.2)]) imply that the same Jacobi matrix corresponds to the operator 
of multiplication by (0— x) / \[Q in the basis {C m /|| C m || } of the space of polynomials. 
It is readily checked that the Charlier weight, viewed as a measure on Z+, satisfies 
the sufficient condition (C), see the proof of Lemma 3.1. Then the same argument 
as in that lemma shows that the space of polynomials in dense in the weighted space 



£ 2 (Z+,Wt 



Charlicr 
9 



) and the above multiplication operator is essentially self-adjoint. 
This is equivalent to saying that the functions C m form an orthonormal basis in 
£ 2 (Z + ) and £)Charhcr j s essentially self-adjoint. Then it follows from the difference 
equation that £)Chariior k as Q m as an eigenvector with eigenvalue The last 

claim of the lemma is now obvious. □ 
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Lemma 4.4. Let s £ R be fixed and assume 6 = 6(N) = N+ (.s + o(l))iV 1 / 2 . Then 
km K%$ff(x, y) = if s Hcrmitc (x, y), x, y £ Z+ . 

Here the assumption on comes from (4.3) and the relation 6 = v/N (Lemma 
4.2). 



Proof. Consider the self-adjoint operators £)C harlior ( wnere 6 > 0) and ZJHcrmite m 
f 2 (Z + ). Let ^o(^+) denote the algebraic linear span of the basis elements S x , x £ 
Z + . By Lemma 3.1 and Lemma 4.3, all these operators have ^o(Z + ) as a common 
essential domain. Moreover, it is evident that if 9 — > oo then £) Charllcr _^ ^Hermite 
on £%{%+). :t follows that Z)Chariicr _^ £)Hormito in the strong reS olvent sense (see 
[RS, Thm. VIII.25]). 

Let us regard K^ llicr and K^ cimitc as operators in £ 2 {Z+). By Lemma 3.1, 

the latter operator is the spectral projection of D Hermlte corresponding to the set 
[s, +oo). By Lemma 4.3, the former operator is the spectral operator of £)Charhcr 
corresponding to the set (4.6). Next, it follows from the description of the spectrum 
of £)Charhor m L emma 4 3 fast instead of the finite set (4.6) we can equally well 
take the continuous interval 

N+l 

+00 



If 6 = 8(N) = N + (s + o(l))Ar 1 / 2 then the left end of this interval can be 
written as s + en where sn — * as N — > 00. Since £) Hermlte has purely continuous 
spectrum, the strong resolvent convergence implies that the spectral projection of 
D'g(N) hcr corresponding to [s+e^r, +00) strongly converges to the spectral projection 

of £) Hermltc corresponding to [s, +00): this is proved exactly as claim (b) in [RS, 
Thm. VIII.24]. □ 

Note that Lemma 4.4 could be obtained from the known asymptotics for the 
Laguerre polynomials [Te] and the well-known connection between the Laguerre 
and Charlier polynomials. 

Lemma 4.4 completes the proof of Theorem 4.1. □ 

5. Another model 

Let N and M be two natural numbers, and (M N ) be the rectangular Young dia- 
gram with N rows and M columns. Given a Young diagram A C (M N ), we denote 
by (M N )/X the skew diagram which is the difference of (M N ) and A. Reading this 
skew diagram from the bottom to the top we get an ordinary Young diagram which 
will be denoted by A: 

(Ai,...,Ajv) = (M — Ajv, . . . ,M — Ai). 

Let tt(m n ) denote the irreducible representation of the symmetric group Snm 
of degree NM, indexed by (M N ). For any n = 0, 1, ... , NM, the restriction of 
7T( M n) to the Young subgroup S n x Snm-u has simple spectrum consisting of the 
irreducible representations of the form n\ ® 7^ (the outer tensor product of the 
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irreducible representations of S n and Snm-u, indexed by A and A, respectively). 
Indeed, this follows from the fact that the skew Schur function S(m«)/a equals the 
ordinary Schur function , as is readily verified using the Jacobi-Trudi formula 
(see [Ma, Ch. I, (5.4)]). 

Let Y n (N, M) stand for the set of Young diagrams with n squares, contained in 
the rectangle (M N ), n = 0, 1, . . . , NM. The above claim shows that the following 
expression defines a probability measure on Y n (N,M), which will be denoted as 

... dim A -dim A 

Mn,N,MW = dMM N } > S YnWM). 

It turns out that if the triple of parameters n, N, M goes to infinity in an appro- 
priate way then the boundary of the random Young diagram distributed according 
to the measure M nj jv,M, after a suitable scaling, tends to a nonrandom curve: This 
is a particular case of the results in Biane [Bil, Thm. 3.1.2] and Pittel and Romik 
[PR, §1.1 and §1.5]. Biane's approach uses free probability. The derivation of Pittel 
and Romik of the explicit form of the limit curve is based on the variational prin- 
ciple. Here we sketch a simple alternative argument. It does not rigorously prove 
the existence of the limit curve but allows one to guess what it is. 

In what follows we will assume M = N and abbreviate M n ^ N = M n ^ N _ N . The 
case of a rectangle (M N ) can be handled in a similar way. We stick to the square 
case M — N to simplify the notation only. 

Assume that N and n go to infinity in such a way that n <~ pN 2 , where p € (0, 1) 
is a fixed parameter. Instead of M n ^ we will be dealing with a modified measure, 
which is obtained by a mixing procedure similar to poissonization: all values n — 
0,1,..., iV 2 are mixed by making use of the binomial distribution on {0,1,..., TV 2 } 
with parameter p. Like the Poisson distribution, the binomial distribution possesses 
the concentration property: as N gets large, the main contribution comes from those 
n's which are close to pN 2 . Thus, one may believe that mixing does not affect the 
asymptotics. 

The resulting measure lives on the set Y(N, N) of all Young diagrams contained 
in (N N ) (no constraints on |A| are imposed): 

M™£(\) = ^pM(l-p) N2 -MM WtN (\), A e Y(N, N). (5.1) 

The next claim, which is similar to Lemma 4.2, shows that the measure M^ff leads 
to the TV-particle Krawtchouk ensemble. 

Denote by ^Krawtchouk j. ne weight function of the Krawtchouk orthogonal poly- 
nomials on the finite set of integers {0,1,..., L} and depending on the parameter 
pe(0,l): 

^Krawtchouk /„\ ( ^ \ nn xf-t \L—x 



^Krawtchouk ^ = ^jp x (l - p) L ~ x , X = 0, 1, . . . , L. 

Lemma 5.1. Under the correspondence A <-» {xi, . . . , xjy} defined by (2.2), random 
Young diagrams A e Y(N, N) distributed according to the measure M^ 1 ^ turn into 
random N -particle configurations in {0,1,..., L} with L = 2N — 1 and such that 

N 

Prob({x 1; . . . , x N }) = const(p, N) ■ J] W^ L Mouk (x t ) ]J (a* - Xj ) 2 . (5.2) 

i=l l<i<j<N 
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Proof. Recall that 

,, . dim A • dim(A) 

Applying Frobenius' dimension formula (4.5) and using (5.1) we obtain the desired 
expression. □ 

The next result describes the limit behavior of the Krawtchouk ensemble "in the 
bulk". 

Fix a real number c such that |c| < 2^Jp{l — p). 

Let apf be an arbitrary sequence of integers such that ~ cN. Given a random 
configuration {xi, . . . , xn} of the iV-particle Krawtchouk ensemble (5.2), we shift it 
by N + aN to the left and obtain in this way a new random iV-particle configuration 
{x' 1; . . . , x' N } C Z: 

x\ = Xi - N - a N , l<i<N. (5.3) 

Proposition 5.2. Under the above assumptions, the random configuration (5.3) 
converges as N — > oo to the translation invariant determinantal random point pro- 
cess on Z with the correlation kernel 

^discrete ^( ^ siny^-y)) z (54) 

?r(x - y) 



where 



( c(l - 2p) \ 
</9 = arccos — -^^^^^^^^= . (5.5) 



e(l-2p) 



2 x /(l-c 2 )p(l-p) 



< 1, so that 



The assumption \c\ < 2y/p{l — p) just means that 
(p is well defined. 

The kernel (5.4), called the discrete sine kernel, first appeared in connection 
with the Plancherel model, see [BOO]. This kernel is a lattice counterpart of the 
celebrated sine kernel 

tt{u — v) 

The result of the proposition is a manifestation of a general phenomenon studied 
in [BKMM]: The discrete sine kernel is the universal correlation kernel arising in 
the bulk limit of discrete orthogonal polynomial ensembles. 

Sketch of proof . The argument is very similar to the proof of Theorem 4.1. We give 
the formal computation below; the justification is omitted. 

Let K m (x;p, L) denote the Krawtchouk polynomial of degree m (necessary in- 
formation about these polynomials can be found in [KS, §1.10]). The normalized 
functions 

K m (x;p,L) = (W^r tch ° u Hx)) 1/2 WKU-^m^K^x-^L) 

form an orthonormal basis in the £ 2 space on the finite set {0,1,..., L}. The 
iV-particle Krawtchouk ensemble (5.2) is a determinantal point process with the 
correlation kernel 

JV-l 

^Krawtchouk {x ^ y) = £ K m (x]p, L)K m (y;p, L). (5.6) 

m=0 
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The Krawtchouk polynomials K m (x;p, L) satisfy the difference equation 
p(L - x)K m (x + l;p, L) + x(2p - l)K m (x;p,L) + x(l - p)K m (x -l;p,L) 

= (pL - m)K m (x; p, L). 
In terms of the normalized functions the difference equation takes the form 

1) 



y/(L- X )(x + l) 



K m (x + l\p,L) + 



x(2p 

Vp(i-p) 



K m (x;p, L) 



y/(L-X+l) X 



K m (x - L) 



pL — m 



-.K m (x;p,L). 



L ' L^pJT-p) 

Let D denote the difference operator defined by the left-hand side of this equa- 
tion. The correlation kernel (5.6) corresponds to the projection on the following 
part of the spectrum of the operator D: 



pL — m 



P) 



0,. 



,N-1 



(5.7) 



{lVpTT 

Recall that L = 2N — 1 and x s=s N + cN + x'. For large N, the three coefficients 
of our difference operator are approximately equal to 



(l + c)(2p-l) 
2 v ' ' ' 2^/pJT^pj ' 

and the set (5.7) approximates the interval 

2p - 1 2p 



Thus, in the limit N 



2a/p(1-p) ' 2Vp(l-p) 
oo we get the difference operator 



(5.8) 



\V^f{x> + 1) + % + J^_-? fV) + \VT=*fV - 1), 



.x' e z, 



and the spectral projection corresponding to the interval (5.8). 

Subtracting the scalar operator / ^ an< ^ ^ v ^ m S 

finally arrive to the difference operator 
f{x' + 1) + /(x' 

and the spectral interval 

c(l - 2p) 



2-y/p(l-p) 

l), a/ez, 



c(l - 2p) + 1 



V(l-c 2 )p(l-p) 



(5.9) 



(5.10) 



V(l-c 2 )p(l-p)_ 

The corresponding spectral projection is given by the discrete sine kernel (5.4). 
Indeed, to study the difference operator (5.9) it is convenient to make the Fourier 
transform from l 2 (Z) to the L? space on the unit circle \z\ = 1. Then (5.9) becomes 
the operator of multiplication by the function z + z = 23?z. Hence, we see that 
our operator has purely continuous (double) spectrum ranging from —2 to 2. It is 
readily seen that the right end of the interval (5.10) is always > 2, while the left end 
is somewhere inside this interval (here we use the assumption \c\ < 2^/p(\ — p) ). 
Thus, in the L 2 space on the circle, our spectral projection becomes the operator of 
multiplication by the characteristic function of the arc going in the counterclockwise 
direction from e~ lip to e z(p , where ip is given by (5.5). In the £ 2 (Z)-realization, this 
is the integral operator with the discrete sine kernel (5.4). □ 
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Corollary 5.3. Let A S Y(AT, iV) 6e i/ie random Young diagram distributed ac- 
cording to the probability measure M^y. Assuming that the boundary of X in the 
scaled coordinates x — i/N , y = j/N has a nonrandom limit described by a curve 
x + y = F(y — x) we can explicitly find F from the equation 

1 - F'(c) ld 1 / c(l - 2p) 
= — = — arccos , 

2 7T 7T \2y/(l-<?)p(l- P ) 



where c ranges over the interval (—2y // p(l — p), 2y / p(l — p) ). 

An additional condition is that the area bounded by the limit curve and the 
coordinate axes in the (x, y) plane has to be equal to p. 

Idea of proof. We observe that the density function (that is, the first correlation 
function) of the point process with discrete sine correlation kernel is the constant 
tp/ft. Then we use the same argument as in item (b) of §1 (see also [BOO, Remark 
1.7]). □ 

One can verify that this result agrees with the formulas in [PR] . The endpoints 
of the interval (—2y/p(l — p), 2y/p(l — p) ) correspond to the endpoints of the 
limit curve that lie on the coordinate axes. 

Note that the result of Proposition 5.2 can be obtained using asymptotics of 
Krawtchouk polynomials obtained in [IS]. The case p = 1/2 is also handled in 
[Jo2, Lemma 2.8]. These papers contain much finer results on the asymptotics but 
obtaining them requires substantially more work. 
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